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We clarify an important aspect of density functional theories, the broadening of the derivative 
discontinuity (DD) in a quantum system, with fluctuating particle number. Our focus is on a 
correlated model system, the single level quantum dot in the regime of the Coulomb blockade. 
We find that the DD-broadening is controlled by the small parameter T/U, where F is the level 
broadening due to contacting and U is a measure of the charging energy. Our analysis suggests, 
that Kondoesque fluctuations have a tendency to increase the DD-broadening, in our model by a 
factor of two. 
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Introduction. Over the years the density functional 
theory (DFT) developed into an important tool to 
study transport properties of nano-systems and sin- 
gle molecules This development occured despite of 
the fact, that often the results are quantitatively sen- 
sitive to the approximations made for the exchange 
correlation (XC) functional, VxcM) underlying such 
calculations.— ~— In principle, dc-transport calculations 
should combine either long-time evolution of wavepackets 
or a KS-bascd quasistationary formalism, in both cases 
with special dynamical XC-functionals.— d 1 : 12 In practice, 
the available ground state functionals are being used. 

The neglect of dynamical correlations for simulation of 
dc-transport was justfied for spinless systems for which 
a Friedel-sum rule holds.— In such systems, approxima- 
tions to the XC-potential of the ground state introduce 
the largest numerical error in the regime of Coulomb 
blockade (CB) where the system ("quantum dot", QD) 
is only weakly coupled to the electronic reservoirs and 
the filling is close to an integer. Its signature is an addi- 
tion energy, U, that largely exceeds the (single particle) 
level spacing of the QD. 

In this Letter we exploit the observation that CB 
is intrinsically an equilibrium phenomenon even though 
it is mostly discussed in its effect on transport 
measurements^; in a broader context it is a typical man- 
ifestation of the reduced compressibility, dn(r)/dfi, of re- 
pulsively interacting fcrmion gases. Therefore, it has a 
reincarnation in XC-functionals of DFT where it appears 
as the derivative discontinuity (DD). 

Starting from the seminal work by Perdew et al.— , the 
DD was almost exclusively discussed in the limiting case 
of decoupled quantum dots, i.e. closed systems. There, 
the XC-functional jumps discontinuously when tuning 
the particle number, N, of a closed system in its ground 
state through an integer value 
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where fj, denotes the electrochemical potential of the 
N particle system which (up to a sign) equals the 
workfunction i 17 : 18 The relation between A and Axe is 
easy to see. Due to Janak's theorem the energy of a 
KS-orbital, index M, in the TV-particle system, ejjj- is 
related directly to the work functions: n N ~ SN =e^~ SN , 
With Eq. ([5J we conclude 
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where Aks is the energy spacing between the lowest un- 
occupied (M = N + 1, LUMO) and the highest occupied 
(M = N, HOMO) KS states, A KS =e%ll N -e%~ 5N , and 
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It follows that the DD, Axe, accounts for the difference 
between the bare, single particle gap, Aks, and the ad- 
dition energy, A, for supplying one more particle. This 
extra energy cost, Axe, related to the repulsive interac- 
tion of fermions confined in a narrow region of space is 
also the origin of the CB and incompressibility. The DD 
in closed iV-particle systems and ways to include it into 
approximate schemes have been a subject of intense re- 
search, recently^ In local or semi-local approximations 
of XC-functionals artifacts in the description of charge 
transfer and transport processes arise, because the DD is 
not accounted for^i - - 2 ^^ 

In open systems the understanding of the DD is still 
relatively poorly developed. In particular, its fate in a 
situation with weakly coupled subsystems, e.g., QD and 
electronic/thermal reservoirs, has not yet been studied 
systematically. This is our focus here. We investigate 
the equilibrium compressibility, dn(r)/dfi, in a generic 
model system, the Anderson (or single site Hubbard) 
model with a repulsive on-site interaction U. The smear- 
ing of the discontinuity ([T} , as a consequence of coupling 
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to a reservoir, can be observed in simulations employing 
the density matrix renormalization group (DMRG) ? . In- 
sight about paramctrical dependency is drawn from an- 
alytical results that also allow the construction of "toy" 
XC-functionals to study CB in DFT. 

Specifically, we report the compressibility dN/d^x near 
integer filling, N — 1, in three different temperature 
regimes. We summarize our findings, (i) For the isolated 
level we have at nonzero temperatures (j3 = 1/T) the ex- 
act result: ^\ N=1 =Te pu / 3 +T-U/2. The derivative 
is finite as T > and the discontinuity is broadened 
by temperature. A dimensionless measure of this smear- 
ing is provided by the number of particles SN that need 
to flow into the level in order to fascilitate the increase 
("jump") of A(N) by A XC « U: 



SN = U dN/d^\ N=1 = /3U(e 0U / 2 +l) 
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In this perspective, the DD is the statement that in the 
zero temperature limit the amount of particles needed 
to drive the jump becomes arbitrarily small, (ii) In the 
presence of a weak coupling to an electronic reservoir, the 
single level acquires a width V . Below a certain cross-over 
temperature, Tp T, we witness that SN stops to decrease 
and the lifetime broadening leads to an intermediate sat- 
uration; we obtain 

SN « 4T/irU(l + 0(T/U)) T K « T < T r (6) 

(iii) At even lower temperatures, below the Kondo scale, 
Ik, where the Abrikosov-Suhl resonance if fully devel- 
oped, we cite an exact asymptotic result: 



SN = 8T/ttU{1 + 0(T/U)) T<T k 
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In cases (ii) and (iii) model XC-functionals arc given, 
that reproduce the CB-features on a qualitative level. 

Anderson model. The Anderson modeled describes a 
single level QD coupled to a reservoir (1Z): 
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where Hn= Eo-k e k<4. k c<7k, e k = — 2tcos(ka) (a: lattice 
spacing) and in the presence of spin-rotational invariance 
(n ff =4d ff ,7V=n t + n;, (h a )=N/2): H QB =e d N + Un t n l 
The observable of interest in DFT is the local density 
and its variation with exernal parameters, e.g. ed and /z: 
N(fj) and dN/dfj, especially near integer fillings, N = 1. 
The inverse, n(N), will then be related to the exchange- 
correlation potential on the QD via 



Vxc(N) = f x(N)-UN/2-e d ; 
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Hartrec and on site potential here must be subtracted. 
To calculate physical observables in the presence of reser- 
voirs, thermal Green's functions provide a convenient for- 
malism. For the thermal occupation numbers we have 
quite generally 



N/2 = (n a ) = T^Q a {\u m ) = J dEf E A a (E) 
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where we have introduced the spectral function A a (E) 
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a. Thermal coupling. In the absense of particle fluc- 
tuations, V = 0, the spectral function can be calculated 
exactly M It is given by (a = —a) 

A a {E) = {\-{n s ))5{E-e d ) + {n s )5(E-{e A + U)) (11) 

The two "Hubbard bands" are reflecting the energy cost, 
U, for adding the second particle to the QD. Recalling 
(JTTJ) we obtain 
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which implies that at integer filling, JV=1, we have 

Hi=£d+U/2. Recalling equation ([9]) we conclude 

^xc|iv=i = 0. By inverting (|12|) we obtain the general 
answer 
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with a(pU,N) = s/i + {eP u ~l)(N-l) 2 . It is implied 
that near integer filling dfi/dN\ N=1 = T + Te 1311 / 2 and 
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From this expression it is obvious how the DD emerges: 
at any nonzero value of the interaction parameter U, the 
slope near N = 1 diverges in the zero temperature limit. 
The divergency occurs in an exponential way because 
those fluctuations in particle numbers that give fi a non- 
vanishing slope are suppressed by a factor of exp/3U/2. 

The diverging slope can also be interpreted in the fol- 
lowing way. At low temperature and near integer filling 
a very small change in the local particle number, SN, 
can increase the effective on-site potential by the finite 
amount U: SN = UdN/dp = /3U(e^ u / 2 + l)" 1 We have 
arrived at Eq. ([5]). 

b. The quantum limit: T — > 0. In the presence of 
the nonvanishing coupling, V > 0, the occupation num- 
bers h a no longer commute with the Hamiltonian, H, and 
the Anderson model becomes nontrivial. Two essential 
changes occur. First, the Hubbard bands acquire a finite 
width, T. As a consequence, the decrease of SN(f3U) 
stops when T falls below T. The residual density fluctu- 
ations near integer filling then are no longer controlled 
by thermal but by quantum fluctuations. The control 
parameter for the latter is T/U; it measures the overlap 
of the Hubbard bands with the Fermi-energy. Second, 
at lowest temperatures, T <c Xk, the spectral function 
acquires a third peak, the Abrikosov-Suhl resonance. 
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FIG. 1: (Color online) The change of the dot filling rescaled 
by | ~ 2 in the Coulomb blockade regime for an on-site repul- 
sion of U — 5t, a hybridization of V — O.St (+) and V — O.lt 
(x) for a M = 30 site system obtained from a ground state 
DMRG. For comparison results for a 50 site calculation are 
also shown (red,blue o) emphasizing convergence with the sys- 
tem size. Data illustrates the scaling of dN/dp with the hy- 



bridization r ~ 
dot occupation 



\V\ . Raw data is shown in the inset. Inset: 



c. Intermediate temperatures: Tk <C T < T,U . We 
first imagine that the Kondo temperature is by far the 
smallest energy scale, in particular Tk "C T, so that the 
Kondo effect can be ignored. This is justified in effec- 
tively finite systems, like large molecules, where the den- 
sity of states (DoS) near the highest occupied molecu- 
lar energy level is not truely continuous. In this situ- 
ation the presence of the reservoirs is easily dealt with 
on a qualitative level by equipping the thermal Green's 
function Eq. (|T0|) with a self energy, E(iu; m ). It deter- 
mines the inverse lifetime T(E) = —Q'S(E) which evalu- 
ates to T(E) = tt\V\ 2 p-ji(E) for noninteracting reservoirs 
(pn(E): reservoir DoS). The qualitative effect can be 
studied in the simplest approximation where we ignore 
the energy dependency of the self energy ("wide band 
limit"). Then, the spectral function for a given Hubbard 
subband takes the Lorcntzian shape 
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As a consequence of the peak broadening, Eq. ([12"]) 
generalizes, 
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with the zero temperature limit 
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where = e d + Eq. (fTf) reveals that at TV = 1 we 
still have p\ = e* d + U/2 at arbitrary T,T values. The 
formulae (|17|19|) combine into a transcendental equation 
for p(N) and by virtue of Eq. ^ into a functional i>xc- 
For calculating the dN/dp we notice, that p enters F e * 
and i^/i-ej-C/ with the opposite sign, implying that the p 
derivative of the denominator of (fT7| vanishes at N = 1 . 
Hence wc derive 
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Eq. © follows via expansion in T/U. 

d. The Kondo limit: T <C Tk- When the tempera- 
ture decreases down to the Kondo scale, T ~ Tk, 



Tk = 



- 7 r| M -e*|| M -e*-C/|/2C/r 
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(21) 

the Abrikosov-Suhl (AS) resonance starts to build up<2i 
(D: conduction electron bandwidth, c « (D 2 /\p~ e d \\p — 
( *d - U\) 1/2 ; c « 0.29 in the wide band limit, D -> oo.) 
When it is fully developed, T -C Tk, its shape is roughly 
Lorcntzian, ^ A s(T : )~(l/7rr)T K /((£:-^) 2 +T K )- 1 , and it 
adds a third resonance to the spectral function 
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with a normalizing coefficient, p=l/(l+Tk/2r). As writ- 
ten, Eq. (|22p has an artificial feature in the sense that the 
peak values at all resonances coincide: l/7rr. A slightly 
more accurate representation incorporates a change in 
the shape of the Hubbard bands in the Kondo regime 
which we here account for by replacing the original width 
of A a with another one, T — > T K , specified below. With 
this caveat we have (t = J dEf E AAs(E) = T K /2T): 
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Eq. ([23| captures qualitative features of N(p) and it 
can be used to construct an LDA for a Kondo-system. 
Eq. (j23j) suggests that the impact of the AS-resonancc 
on the compressibility is small as Tk/r. The main im- 
pact of Kondocsquc fluctuations comes here from the 
rcnormalization of the shape of the Hubbard peaks. 
Indeed, the exact compressibility known from Bcthc- 
Ansatz calculations^ reads in the limit of large U: 
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implying Eq. ([7]). Analogous to Eq. ©, an estimate 
based on (|23[) would give 4T k /ttU 2 which suggests T K w 
2r when comparing with (|2"4")l . Wc conclude that in the 
Kondo- regime dN/ dp should be enhanced by a factor ~ 2 
reflecting a stronger tendency for charge-fluctuations. 
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e. DMRG- calculation. To illustrate our analytical 
arguments DMRG-calculations have been performed on 
systems including M = 30, 50 sites. The available sys- 
tem sizes do not allow us to resolve the Kondo-scale, yet, 
but the expected finite slope of dN/d/j, at integer filling 
is clearly visible. The data also shows the collaps on a 
single curve when rescaled by |V^|~ 2 . 

Discussion. Our survey of analytical results obtained 
in the symmetric Anderson model suggests that the DD 
is broadened in systems coupled to a reservoir: the par- 
ticle transfer, SN, needed to shift the local XC-potential 
by the on-site interaction energy ~ U is not infinitesi- 
mally small. Within the model considered, the particle 
transfer is a two-parameter function, SN(j3U,T/U), that 
smoothly interpolates between a high temperature and a 
low-temperature ("quantum") regime. 

This result has implications for model studies of quan- 
tum transport within the framework of time depen- 
dent DFT. The importance of the DD for such trans- 
port simulations has been emphasized in several recent 
works . 7 ' 22 ! 25 Our work implies, that effective functionals 
used in such simulations should exhibit a parametric de- 



pendency on T/U. In particular, only terms quadratic 
in the coupling, V 2 , appear in single channel quantum 
transport. The ABALDA- functional proposed in Ref. 
[25| does not adhere to this principle since it depends 
explicitly on the parameter [//Vimk- 

Our results also have implications for molecules, i.e. 
systems where one subsystem couples weakly to a small 
number of other subsystems, but not to a (macroscopic) 
reservoir, proper. In this situation, the spectral function 
Act, Eq. ([TT]) . will translate into the local spectral func- 
tion of the given subsystem: each Hubbard peak acquires 
a splitting indicative of the hybridization of states with 
the environment. Again, the amount of charge needed 
to fill into the subsystem to drive vxc up by U will not 
be zero but rather reflect this hybridization induced sub- 
structure. 
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* Electronic address: ferdinand.evers@kit.edu 

1 S. H. Ke, W. Yang, and H. U. Baranger, Nano Letters 8, 
3257 (2008). 

2 A. Gagliardi, G. C. Solomon, A. Pecchia, T. Frauenheim, 
A. Di Carlo, N. S. Hush, and J. R. Reimers, Phys. Rev. B 
75 (2007). 

3 T. Markussen, R. Stadler, and K. S. Thygesen, Nano Lett. 
10, 4260 (2010). 

4 M. Dell' Angela, G. Kladnik, A. Cossaro, A. Verdini, 
M. Kamenetska, I. Tamblyn, S. Y. Quek, J. B. Neaton, 

D. Cvetko, A. Morgante, et al., Nano Lett. 10, 2470 (2010). 

5 A. Mishchenko, D. Vonlanthen, V. Meded, M. Buerkle, 
C. Li, I. V. Pobelov, A. Bagrets, J. K. Viljas, F. Pauly, 
F. Evers, et al., Nano Lett. 10, 156 (2010) 

6 S. Ke, H. U. Baranger, and W. T. Yang, J. Chem. Phys. 
126, 201102 (2007). 

' C. Toher, A. Filippetti, S. Sanvito, and K. Burke, Phys. 
Rev. Lett. 95, 146402 (2005). 

8 J. R. Reimers, Z.-L. Cai, A. Bilic, and N. S. Hush, Ann. 
N.Y. Acad. Sci. 1006, 235 (2003). 

9 F. Evers, F. Weigend, and M. Koentopp, Phys. Rev. B 69, 
235411 (2004). 

N. Sai, M. Zwolak, G. Vignale, and M. D. Ventra, Physical 
Review Letters 94, 186810 (2005). 

1 G. Stefanucci and C.-O. Almbladh, Phys. Rev. B 69, 
195318 (2004). 

2 S. Kurth, G. Stefanucci, C.-O. Almbladh, A. Rubio, and 

E. K. U. Gross, Phys. Rev. B 72, 035308 (2005). 

3 P. Schmitteckert and F. Evers, Physical Review Letters 
100, 086401 (2008). 

4 H. Bruus and K. Flensberg, Many Body Quantum The- 
ory in Condensed Matter Physics (Oxford Graduate Texts, 
2004). 

5 L. I. Glazman and M. Pustilnik, arXiv:cond-mat/0501007, 



Lecture notes of the Les Houches Summer School 2004 
"Nanophysics: Coherence and Transport", eds. H. Bouch- 
iat et al. (Elsevier, 2005) p. 427 (2005). 
J. Perdew, R. Parr, M. Levy, and J. Balduz, Phys. Rev. 
Lett. 49, 1691 (1982). 

F. P. Roselli, A. B. F. da Silva, and K. Capelle, 
arXiv:physics/061180v2 (2006). 

P. Mori-Sanchez, A. J. Cohen, and W. Yang, Phys. Rev. 

Lett. 100, 146401 (2008). 

J. F. Janak, Phys. Rev. B 18, 7165 (1978). 

C.-O. Almbladh and U. v. Barth, Phys. Rev. B 31, 3231 

(1985). 

S. Kiimmel and L. Kronik, Rev. Mod. Phys. 80, 3 (2008). 
M. Koentopp, K. Burke, and F. Evers, Physical Review 
B (Condensed Matter and Materials Physics) 73, 121403 
(2006). 

F. Evers and K. Burke, arXiv:cond-mat/0610413; Chapt. 
24 in CRC Handbook on Molecular Electronics, ed. E. Ly- 
shevski (2006). 

A. C. Hewson, The Kondo Problem to Heavy Fermions 

(Cambridge Studies in Magnetism, 1995). 

S. Kurth, G. Stefanucci, E. Khosravi, C. Verdozzi, and 

E. K. U. Gross, Phys. Rev. Lett. 104, 236801 (2010). 

C. Fiolhais, F. Nogueira, and M. Marques, eds., A Primer 

in Density Functional Theory (Springer, 2003). 

J. P. Bergfield, Z-F. Liu, K. Burke and C. A. Stafford, 

arXiv:1106.3104. 

Most current density functionals rely upon local (LDA) 
and semi-local (GGA) approximations 2 ^; they exhibit 
a continuous evolution with N which implies that the 
LUMO JV+,5JV and the LUMO^- 5 ^ see essentially the same 
effective potential; therefore A^ A ' GGA = 0. 



